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Abstract 



We have used the first 2600 nontrivial zeros 7; of the Riemann zeta 
function calculated with 1000 digits accuracy and developed them into 
the continued fractions. We calculated the geometrical means of the 
denominators of these continued fractions and for all cases we get values 
close to the Khinchin's constant, what suggests that 7; are irrational. 
Next we have calculated the n-th square roots of the denominators q n of 
the convergents of the continued fractions obtaining values close to the 
Khinchin-Levy constant, again supporting the common believe that 7/ 
are irrational. 



1 Introduction 

Bernhard G.F. Riemann has shown [TB] that the l.h.s. of the identity valid only for 



can be analytically continued to the whole complex plane devoid of s = 1 by means 
of the following contour integral: 



3?[s] > 1: 




s = 



a + it 



(1) 




(2) 



where the integration is performed along the path V 
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+■ 



Now dozens of integrals and series representing the ((s) function are known, for 
collection of such formulas see for example the entry Riemann Zeta Function in [20J 
and references cited therein. 

The C,{s) function has trivial zeros —2, —4, —6, . . . and infinity of nontrivial com- 
plex zeros pi = (3i + iji in the critical strip: (3i G (0, 1). The Riemann Hypothesis 
(RH) asserts that fii — | for all I — i.e. all zero lie on the critical line = |. 
Presently it is added that these nontrivial zeros are simple: C'(Pz) 7^ — many 
explicit formulas of number theory contain ('{pi) i n the denominators. In 1914 G. 
Hardy [6] proved that infinitely many zeros of ((s) lie on the critical line. A. Selberg 
[j~8] in 1942 has shown that at least a (small) positive proportion of the zeros of ((s) 
lie on the critical line. The first quantitative result was obtained by N. Levinson in 
1974 [H] who showed that at least one-third of the zeros lie on the critical line. In 
1989 B. Conrey [2] improved this to two-fifths and quite recently with collaborators 
[lj to over 41%. It was checked computationally [5J that the 10 13 first zeros of the 
Riemann Zeta function fulfill the condition fti — |. A. Odlyzko checked that RH is 
true in different intervals around 10 20 [10], 10 21 [UJ, 10 22 [H], see also [5] for the 
two billion zeros from the zero 10 24 . 

There is no hope to obtain the analytical formulas for the imaginary parts 7; of 
the nontrivial zeros of ((s) but the common belief is that they are irrational and 
perhaps even transcendental [13]. The problem of any linear relations between 7/ 
with integral coefficients appeared for the first time in the paper of A.E. Ingham [7] 
in connection with the Mertens conjecture. This conjecture specifies the growth of 
the function M(x) defined by 




(3) 



where /i(n) is the Mobius function 




for n = 1 

when p 2 \n 

when n = "pvpi . . . p r 



(4) 



The Mertens conjecture claims that 



\M(x)\ < 



(5) 
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From this inequality the RH would follow. A. E. Ingham in [7] showed that the 
validity of the Merten's conjecture requires that the imaginary parts of the nontrivial 
zeros should fulfill the relations of the form: 

N 

£> 7l = 0, (6) 
1=1 

where q are integers not all equal to zero. This result raised the doubts in the 
inequality ^ and indeed in 1985 A. Odlyzko and H. te Riele [15] disproved the 
Merten's conjecture. 

In this paper we are going to exploit two facts about the continued fractions: the 
existence of the Khinchin constant and Khinchin-Levy constant, see e.g. [H §1.8], 
to support the irrationality of 7;. Let 

1 

r = [a (r); ai(r), a 2 (r), a 3 (r), . . .] = a (r) H (7) 

ai(r) H 



02 (r) H 

a 3 (r) + '•• 

be the continued fraction expansion of the real number r, where ao(r) is an integer 
and all a^(r) with /c > 1 are positive integers. Khinchin has proved [8], see also |17j . 
that 

lim ( ai (r)...a n (r))™ = TT \l + — f = K « 2.685452001 .. . (8) 

is a constant for almost all real r jH §1.8]. The exceptions are of the Lebesgue 
measure zero and include rational numbers, quadratic irrationals and some irrational 
numbers too, like for example the Euler constant e = 2.7182818285... for which 
the limit (p]) is infinity. The constant K is called the Khinchin constant. If the 
quantities 

K(r; n) = {a l {r)a 2 {r) . . . a n {r)) n (9) 

for a given number r are close to Kq we can regard it as an indication that r is 
irrational. 

Let the rational p n /q n be the n-th partial convergent of the continued fraction: 

— = [a ;ai,a 2 ,a 3 ,...,a n ]. (10) 

Qn 

For almost all real numbers r the denominators of the finite continued fraction 
approximations fulfill: 

lim (q n (r)) 1/n = e " 2/121n2 = L = 3.275822918721811 . . . (11) 

where Lq is called the Khinchin — Levy's constant [4, §1.8]. Again the set of excep- 
tions to the above limit is of the Lebesgue measure zero and it includes rational 
numbers, quadratic irrational etc. 
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2 The computer experiments 



First 100 zeros 7; of ((s) accurate to over 1000 decimal places we have taken from 
[12] . Next 2500 zeros of ((s) with precision of 1000 digits were calculated using 
the built in Mathematica v.7 procedure ZetaZero[m]. We have checked using 
PARI/GP [12] that these zeros were accurate within at least 996 places in the 
sense that in the worst case \({pi)\ < 10~ 996 , I = 1, 2, . . . , 2600. PARI has built 
in function contf rac(r, {nmax}) which creates the row vector a(r) whose com- 
ponents are the denominators a n (r) of the continued fraction expansion of r, i.e. 
a = [a (r); ai(r), . . . , a n (r)] means that 

r Ki a (r) H (12) 



1 



w) + 



1 



a n (r) 



The parameter nmax limits the number of terms a nmax (r); if it is omitted the 
expansion stops with a declared precision of representation of real numbers at the 
last significant partial quotient. 

By trials we have determined that the precision set to \p 2200 is sufficient in the 
sense that scripts with larger precision generated exactly the same results: the rows 
a(7;) obtained with accuracy 2200 digits were the same for all / as those obtained for 
accuracy 2600 and the continued fractions accuracy set to 2100 digits gave different 
denominators a n {^/i) With the precision set to 2200 digits we have developed the 
1000 digits values of each 7/, I = 1, 2, . . . 2600, into the continued fractions 

7z = [a (l);a 1 (l),a 2 (l),a 3 (l),. . . ,a n( i)(l)] = a(Z) (13) 

without specifying the parameter nmax, thus the length of the vector a(Z) depended 
on 7; and it turns out that the number of denominators was contained between 
1788 and 2072. The value of the product a\a 2 ■ ■ ■ a n (/) was typically of the order 
10 800 — 10 870 . Next for each I we have calculated the geometrical means: 

(nil) \ W 

Ki{n{l))= \ T[a k (l)\ . (14) 

The results are presented in the Fig.l. Values of Ki{n{l)) are scattered around 
the red line representing K$. To gain some insight into the rate of convergence of 
Ki(n(l)) we have plotted in the Fig. 2 the number of sign changes S(l) of Ki(m) — Kq 
for each / when m = 100, 101, . . . n(l), i.e. 

Sk(1) = number of such m that {K t (m + 1) - K )(Ki(m) - K ) < 0. (15) 
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The largest Sk{1) was 122 and it occurred for the zero 7194 and for 381 zeros there 
were no sign changes at all. In the Fig. 3 we present plots of Ki(m) as a function 
of m for a few zeros ji. 

Let the rational p n (i)(ll) / Qn(r\{ll) be the n-th partial convergent of the continued 
fractions (13): 

Pn(i){ll) n , na , 
w = a(Z). 16 

Qn(l)hl) 

For each zero we have calculated the partial convergents p n {i){li) I <in(i){.li) ■ Next 
from these denominators q n (i)(li) we have calculated the quantities Li{n{l))\ 

Li(n(l)) = (q n (i)) 1/n(l) , Z = l,2,...,2600 (17) 

The obtained values of L/(n(/)) are presented in the Fig. 4. These values scatter 
around the red line representing the Khinchin — Levy's constant Lq. As in the case 
of Ki(m) the Fig.5 presents the number of sign changes of the difference Li(m) — L 
as a function of the index m of the denominator of the m-th convergent p m /Qm 

Sl(1) = number of such m that (Li(m + 1) — L )(L;(m) — L ) < 0. (18) 

The maximal number of sign changes was 136 and appeared for the zero 71389 and 
there were 396 zeros without sign changes. 

In the Fig. 6 we have plotted the "running" absolute difference between Ki(m) 
and Kq averaged over all 2600 zeros: 

. 2600 

A «( m ) = 7^Y,\Ki(™) ~ K \, m = 100, 101,... 1788 (19) 
zoOO z — ' 

1=1 

and the similar average for the difference between Li(m) and L : 

. 2600 

Mm) = —J^^im) - L \, m = 100,101,... 1788. (20) 
zoOO z — 4 

1=1 

These two averages very rapidly tend to zero. Although it does not prove nothing, 
the fact that the curves representing Ak{itl) and A^m) almost coincide is very 
convincing. In the inset the plot on double logarithmic scale reveals that both Aiim) 
and Ak(m) decrease like Ck,l/ 'y/rn where Ck = 2.3868 . . . and Cl = 2.4473 .... It 
is a pure speculation linking the power of these dependencies m" 1 / 2 to the ordinate 
of the critical line. 



3 Concluding remarks 

There are generalizations of above quantities K{n) given by ([9]). It can be shown 
that the following s-mean values of the denominators ctfc(r) of the continued fraction 
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for a real number r: 

l/s 



M{n,s;r) = (-£(a ft (r))') (21) 



k=l 



are divergent for s > 1 and convergent for s < 1 for almost all real r [H §1.8]. It 
can be shown that for s < 1 



lim M(„, s; r) = g log 2 (l - |^_) J 



1/s 

= K s (22) 



where M(n, s; r) for almost all r are the same. The quantities (21 ) can be computed 
for imaginary parts of nontrivial zeta zeros M(n, s; 7;) and compared with values of 
K s but we leave it for further investigation. 

The continued fraction were used in the past in the Apery's proof of irrationality 
of C(3). In the paper [3] values of ((n) for all n > 2 were expressed in terms of 
rapidly converging continued fractions. These results were analytical, but in case 
of the nontrivial zeros of the ((s) function we are left only with the computer 
experiments. The results reported in this paper suggests that they are irrational. 
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Fig.l The plot of K^nil)). 
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Fig. 2 The number of such m that (Ki(m + 1) — K )(Ki(m) — K ) < for each I 
(the initial transient values of m were skipped— - sign changes were detected for 

m = 100, 101,... n(Z)). 
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Fig.3 For 7 166 there was no sign change of the difference iT 166 (m) — K . For 7 194 
there were 122 sign changes of the difference Ki 9i (m) — K — it was the largest 
number of sign changes among all zeros. For 71434 there were 63 sign changes of 

the difference ^1434 (m) — K . 
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Fig.4 The plot of Lz(ra(Z)). 
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;.5 The number of such m that (Li(m + 1) — L )(Lz(m) — L ) < for each I (the 
initial transient values of m were skipped — sign changes were detected for 

m = 100, 101,... n(Z)). 
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Fig. 6 The averaged over all 2600 zeros differences \Ki(m) — K \ and \Li(m) — L \ 
plotted for m = 100, 101, . . . 1788). In the inset the same curves are plotted on the 
double logarithmic scale together with fits obtained by the least square method. 
The equations of the fits are 2.3868/m ' 5019 for A^im) (green line) and 
2. 4473 /m - 5028 for Al (to) (blue line) 
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